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INTRODUCTION 


Underlying  the  measurement  of  the  upwelling  radiation  from  a ground  or  sea 
surface  or  a cloud  layer  by  a meteorological  satellite  such  as  SMS/GOES  1 
is  the  theory  of  atmospheric  radiative  transfer.  Without  a reasonable 
understanding  of  the  ramifications  attendant  to  this  theory,  interpretation 
and  analysis  of  the  radiance  data  from  the  satellite,  either  in  streams  or 
imagery,  would  be  inadequate  or  completely  erroneous.  The  importance  of 
the  SMS/GOES  series  to  many  Army  missions  such  as  air  operations,  target 
acquisition,  terminal  homing,  and  the  detection  and  tracking  of  severe 
weather  was  already  noted  in  the  first  report,  "Mesoscale  Determination  of 
Cloud-Top  Height:  Problems  and  Solutions,"  and  will  not  be  reiterated. 

The  present  report,  the  second  of  a proposed  series  of  papers  on  the  theo- 
retical and  experimental  applications  of  satellite  observations,  will  deal 
in  a more  understandable  manner  with  the  theory  of  radiative  transfer  as 
applied  to  a cloudy  atmosphere.  Next,  the  report  will  show  how  the  equation 
of  radiative  transfer  may  be  solved  for  a nonhomogeneous  atmosphere,  with 
brief  mention  of  some  other  techniques.  Since  this  equation  is  to  be  applied 
to  cloudy  atmospheres,  a brief  exposition  of  cloud  microstructure  is  then 
in  order  with  particular  attention  to  stratiform  and  cumuliform  clouds. 

Finally,  a detailed  examination  will  be  made  of  under  what  microphysical 
conditions  a cloud  may  be  unmistakeably  regarded  as  a blackbody  radiator 
and  what  this  implies  in  the  determination  of  cloud-top  heights  and  in 
retrieval  of  surface  temperatures. 


THEORY  OF  RADIATIVE  TRANSFER 

There  are  several  well-known  books  on  radiative  transfer  [1,2,3]  and  a vast 
amount  of  open  literature  and  technical  reports  too  numerous  to  cite  here. 
Chandrasekhar's  book  [1]  is  prepared  primarily  for  the  astrophysicists  in 
their  investigation  of  stellar  atmospheres;  Goody  [2]  provides  excellent 
reference  materials  on  the  calculation  of  atmospheric  transmission  functions; 
Kondratyev  [3]  is  more  comprehensive  and  will  be  found  useful  in  meteorological 
applications,  though  less  rigorous  in  his  treatment  of  radiative  transfer. 

In  this  short  section,  the  author  summarizes  the  experiences  gained  in  his 
study  of  radiative  transfer  and  hopes  to  present  the  theory  in  a lucid 
manner  so  that  anyone  who  has  been  barely  exposed  to  the  problems  of  radia- 
tive transfer  may  be  able  to  follow  the  exposition.  There  will  be  nothing 
new  or  fresh  in  this  approach  to  the  theory;  many  of  the  ideas  are  drawn 
from  previous  literature. 

Radiative  transfer  theory  is  a quantitative  study  of  the  transfer  of  radiant 
energy  through  a medium  which  can  scatter,  absorb,  and  emit  radiation.  In 


representing  radiative  energy,  there  are  two  often-used  radiometric  quantities: 
radiance  and  irradiance  (or  emittance).  In  astrophysics,  they  are,  respectively, 
specific  intensity  and  specific  flux.  Radiance  or  intensity  is  defined  as  the 
energy  radiated  from  a source  per  unit  time  per  unit  projected  area  per  unit 
solid  angle  per  unit  spectral  interval.  In  satellite  meteorology,  when  referred 
to  the  infrared  band,  the  physical  unit  radiance  is  often  given  in  milliwatt 
per  square  meter  per  steradian  per  unit  wavenumber.  Spectral  irradiance  (or 
emittance)  or  flux  is  defined  as  the  energy  irradiated  upon  or  emitted  from  a 
unit  surface  per  unit  time  per  unit  spectral  interval.  The  physical  unit  of 
irradiance  is  milliwatt  per  square  meter  per  unit  wavenumber.  In  the  solar 
radiation  region,  it  is  usually  expressed  in  calorie  per  minute  per  square 
centimeter  per  unit  wavelength.  The  total  flux  density  of  solar  radiation, 
which  is  given  in  calorie  per  minute  per  square  centimeter  or  langley  per  min- 
ute, as  used  in  meteorology,  is  obtained  by  integrating  solar  irradiance  over 
all  wavelengths. 

With  these  definitions,  the  equation  of  radiative  transfer  along  an  arbitrary 
path  according  to  the  law  of  energy  conservation  can  be  written 


dl. 


■ -Vx  ds  * h 


ds 


(1) 


where  I denotes  the  monochromatic  radiance, a function  of  position,  direction, 

A 

and  wavelength  X;  Bg  the  volume  extinction  coefficient,  a function  of  wavelength; 

s is  the  path  length;  and  the  source  function  coefficient  to  be  determined, 

also  a function  of  position,  direction,  and  wavelength.  The  first  term  on  the 
right  side  of  Eq.  (1)  represents  the  reduction  of  radiation  energy  due  to  the 
extinction  process  in  the  medium,  whereas  the  second  term  denotes  the  increase 
of  radiation  energy  carried  by  scattering  and  emission  processes. 

The  volume  extinction  coefficient  8g  (cm-1)  is  related  to  the  mass  extinction 
coefficient  kg  (g_1  cm2)  by 


Be  = p ke 


(2) 


where  p is  the  density  of  the  medium  in  which  the  extinction  of  radiation 
takes  place.  In  general,  the  extinction  coefficient  consists  of  two  radia- 
tion components;  namely,  the  volume  (or  mass)  scattering  6 and  absorption 
coefficient  8,.  Thus, 


6e  = Bs  + Ba 


(3) 
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Normally,  molecules  and  particulates  in  the  atmosphere  scatter  and  absorb 
solar  and  infrared  radiation.  In  the  infrared  window  of  9 to  12  micrometers 
which  is  of  interest  in  this  report,  the  author  is  not  too  concerned  with 
the  scattering  by  gas  molecules.  In  comparison  with  the  magnitude  of  the 
scattering  by  aerosol  and  cloud  particles,  molecular  scattering  is  entirely 
negligible  in  the  infrared  spectrum.  In  the  solar  spectrum,  however,  molecular 
and  particulate  scatterings  are  both  important. 

Returning  to  Eq.  (1),  note  that  in  the  absence  of  an  external  source  the 
equation  reduces  to  a form  which,  upon  integration,  gives  one  of  the  funda- 
mental laws  in  radiation,  that  is,  Beer's  or  Bouguer's  law  of  transmission. 

The  geometry  is  now  fixed  to  avoid  an  arbitrary  path  length.  Figures  1 and  2 
depict  such  a geometry  in  two  and  three  dimensions.  The  radiance  becomes  a 
function  of  optical  path  or  depth  x,  zenith  angle  e,  and  azimuth  angle  <f>. 
Moreover,  the  atmosphere  has  inadvertently  been  treated  as  being  made  up  of 
many  stratified  layers,  that  is,  a plane-parallel  atmosphere  with  infinite 
horizontal  extent.  As  a matter  of  fact,  this  is  how  the  planetary  atmosphere 
is  being  handled  in  the  vast  majority  of  publications  on  visible  and  infrared 
radiative  transfer  either  in  a clear  or  a cloudy  medium.  From  the  figure, 

Eq.  (1)  can  be  transformed  into  the  following  expression: 


dlx  (t  ; e,  <f>)  = -B  Ix(x;  e,  <f>)dz/cos  e + jx(x;  $ , <j>)  dz/cos  6 , (4) 

from  which  we  arrive  at  a form  of  the  radiative  transfer  equation  usually 
found  in  the  literature,  by  setting 


dx  = -3edz  , (5) 

and  y = cos  e;  thus, 


dl  (x;  e,  <f>) 

v — : dx = 05  ^ 05  ^ * (6) 

wnere  the  source  function  J,  = j,/8„. 

As  it  stands,  the  above  differential  equation  can  be  solved  analytically 
without  having  to  know  what  the  source  J stands  for.  In  the  infrared  region, 
a nonscattering,  plane  parallel  atmosphere  which  is  in  local  thermodynamic 
equilibrium  may  be  considered.  Thus,  the  source  in  this  case  is  simply  given 
by 


5 


(7) 


0 A ( t ; 6,  <j>)  = Ba(T)  , 


where  B (T)  is  the  well-known  Planck  function.  In  the  wavelength  domain, 

A 


it  is  given  by 


B , ( T ) = 


whereas  in  the  wavenumber  domain 


^5(ehc/ktA_1 j ’ 


p m _ 2hv3c2  , 

vTj  ■ ; hcv/kT' ; ’ 1 

(e  -1) 

where  h is  the  Planck  constant,  c the  speed  of  light,  k the  Boltzmann  con- 
stant, T the  temperature  of  the  medium,  and  v the  wavenumber  (cm-1). 
According  to  Rossini  [4],  these  constants  are,  in  cgs  units, 

2hc2  = 1.191062  x 1011  ergs  sec-1  cm-/  sr-1  pm4, 
hc/k  = 1.438786  x 104  pm  °K, 

c = 2.99792458  x 1010  cm  sec-1  (in  vacuum).  ( 


For  a purely  emissive  atmosphere  in  local  thermodynamic  equilibrium,  the 
radiative  transfer  Eq.  (6)  becomes 


dl x(t  ; 6 , <j>) 


= I (t  ; 6 , <f>)  - B,  (T)  , 


Now  the  contribution  to  extinction  comes  solely  from  absorption.  This  equation 
holds  for  blackbody  radiation  and  can  be  readily  solved  by  introducing  the 

integrating  factor  e-t//>J.  In  the  retrieval  of  temperature  profiles  from  the 
satellite  vertical  temperature  profile  radiometer  (VTPR)  data,  Eq.  (11)  is 
solved  in  a reverse  manner  to  determine  the  Planck  function  from  which  the 
atmospheric  temperature  is  extracted.  In  theory,  this  equation  governs  the 
generation  of  infrared  imagery  in  clear  air. 

The  additional  contribution  to  the  source  function  due  to  multiple  scattering 
is  found  for  a scattering  atmosphere. 


, <t>'  )sin  e1  de 1 d4> 1 . 


y*  2tT 

Ja(t;  9»  <i>)  = 47  J J pa(t;  e,  4>;  e',  $')  ix(T;  e' 

(12) 

It  is  this  source  function  that  complicates  the  solution  of  the  radiative 
transfer  equation.  In  the  above  equation,  p (x;  0,  <j>;  6',  ^ ' ) is  the 

A 

normalized  phase  function,  representing  the  fraction  of  incident  radiant 
energy  scattered  by  a unit  volume  at  a level  t in  the  direction  of  (e,  4>) 
when  that  volume  is  illuminated  from  the  direction  (o',  <j>').  In  representing 
the  scattering  of  polarized  light,  the  intensity  components  can  be  described 
by  the  four  Stoke 's  parameters.  The  phase  function  now  consists  of  four-by-four 
elements  and  is  called  the  phase  matrix.  The  Stoke's  parameters  make  it 
possible  to  determine  the  state  of  polarization,  which  is  of  great  value  in 
explaining  many  atmospheric  optical  phenomena  but  of  little  concern  in  the 
infrared  spectral  region  where  emission  dominates. 

To  simplify  Eq.  (12),  use  u(=  cos  0)  as  the  functional  variable.  Equation 
(6)  then  assumes  the  following  form: 


dl  (t  ; n , 4>) 


= I x (t  ; v , 4> ) 


r2trr+] 

/, 


P,(t;  u,  u1 , *')  I(  ; 1 , *' ) d/  d>‘  ; (13) 


where  Z = 
path.  0 


yBe  is 


the  albedo  for  single  scattering,  a function  of  optical 


Equation  (13)  is  valid  for  solar  radiation  in  the  visible  band.  However, 
the  sun's  rays  incident  upon  the  plane-parallel  atmosphere  are  generally 
taken  to  be  coming  in  a parallel  beam  and,  hence,  are  azimuthal -dependent. 
Solar  radiation  is  different  from  thermal  radiation  in  that  the  latter  may 
be  considered  isotropic,  and,  hence,  azimuthal-independent.  A parallel 
beam  of  solar  radiation  incident  upon  the  atmosphere  or  the  earth  in  the 
direction  (-yo,<j>o)  has  an  angular  distribution  expressible  in  terms  of 

Dirac's  delta  function  in  the  form 


I . (x  : 
mcv 


*)  = 


Es  (x ) 


6(y  - 


"o> 


<5(<t>  - O 


(14) 


n.  v n*- 


It  is  convenient  to  separate  the  direct  solar  radiation  field  from  the 
diffuse  radiation  field  in  this  analysis.  Thus,  the  total  spectral 
radiance  or  intensity  may  be  written  as 


1a(t;  y,  $)  = Iinc(T;  y,  <j>)  + I di  f ( t ; y,  <t>) 


(15) 


which  may  then  be  substituted  into  Eq.  (13)  to  yield 
u — dif(^l  M.»  ll  = Idif(T;  y > <f> ) - p(x ; y,  <f>;  -y^,  4>p)Es(x) 


4tt 


p(r  ; y , 


p'  > <f>‘  ) Idi1:(T  ; y'  , 4) ' ) dy  'd<J>'  , 

(16) 


where  E$(t)  = TTFe'r^tJo.  In  this  equation,  the  wavelength  subscript  A has 

been  dropped  since  it  is  felt  that  no  misunderstanding  should  arise  by  now 
as  to  the  monochromatic  nature  of  the  radiative  transfer  equation  as 
presently  formulated.  The  subscript  dif  will  also  be  dropped  since  only 
the  diffuse  radiation  field  will  be  dealt  with.  However,  the  solar  term 
Es,  which  is  known  as  the  reduced  incident  or  direct  radiation  in  the 

literature,  must  be  included  in  the  final  accounting  of  reflected  and 
transmitted  radiation. 


In  theory,  Eq.  (16)  is  the  governing  equation  for  visible  imagery.  It 
applies  to  a scattering  and  absorptive  atmosphere,  whereas  Eq.  (11)  applies 
to  an  absorptive  and  emissive  one.  An  atmosphere  which  absorbs  necessarily 
emits.  The  emission  occurs  at  much  longer  wavelength  or  small  wavenumber, 
the  so-called  long-wave  radiation  in  meteorology.  In  a purely  scattering 
medium,  = 1 . In  a purely  absorptive  and  emissive  one,  = 0.  In  general, 

in  an  atmopshere  where  scattering,  absorption,  and  emission  take  place 
simultaneously,  Eq.  (11)  is  added  to  Eq.  (16)  to  obtain 


dI(M  u.  1 ) = 


I (i '»  Pi  4>)  - (1  - mQ)  Ba[T(t)] 


“T/Un 

4-  Fe  P ( t ; y , <f> » - Uq , <tg) 


5Q  /-2Y  + 1 

4 ~ / I p(t;  y,  <t>;  y1,  <}>')  I(x;  y',  <j>' ) dy ' d<f> 1 
Jo  */-l 


(17) 
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This  radiative  transfer  equation,  in  its  full  ominous  appearance,  describes 
the  monochromatic  transfer  of  solar  radiation  in  a vertically  inhomogeneous 
atmosphere  which  is  assumed  to  be  in  local  thermodynamic  equilibrium.  This 
integro-differential  equation  is  difficult  to  solve  for  the  general  case  of 
multiple  scattering;  thus,  numerical  te<.nniques  are  normally  being  employed 
Since  solar  and  thermal  infrared  spectra  are  in  two  distinct  regions  where 
overlapping  can  be  ignored,  solar  and  in'rared  radiation  can  be  treated 
separately  [5];  thus,  in  the  infrared  re  on  the  third  term  on  the  right 
side  of  Eq.  (17)  can  be  neglected  and  t following  is  left 

)] 


')  I ( x ; u',  *')  du'd*'  (18) 


which  is  what  is  needed  in  the  determination  of  the  radiation  properties 
of  several  cloud  types.  Owing  to  the  isotropic  radiation  pattern  emitted 
by  the  earth's  surface,  we  may  assume  to  a good  approximation  that  the 
emergent  radiation  is  independent  of  azimuth  angles,  thereby  greatly 
improving  the  appearance  of  the  phase  function  and  hence  the  transfer 
equation.  Moreover,  if  the  choice  is  to  work  with  relatively  thin  layers 
so  that  each  layer  becomes  essentially  homogeneous  in  its  microstructure, 
then  the  single  scattering  albedo  and  the  phase  function  are  no  longer 
functions  of  the  optical  path  or,  more  appropriately,  optical  depth  in 
view  of  the  plane-parallel  medium.  In  doing  so,  Eq.  (18)  is  reduced  to 


u ^hir1=  Ur;  p)  - 0 - SQ)  Bv[T(t)] 


p(y.  u' ) 


I ( t ; u')  du' 


where  the  phase  function  is  given  by 


p(p. 


p(u,  <f>;  w'  > <f>' ) d<t>' 


(19) 


(20) 


This  is  the  form  of  the  radiative  transfer  equation  usually  given  in  the 
literature  for  infrared  radiation.  The  variable  t is  also  used  to  denote 
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the  levels  at  which  upward  and  downward  radiation  emerges.  The  emergent 
radiances  thus  illuminate  two  hemispheres;  therefore,  it  would  be  convenient 
to  separate  the  two  hemispheres  so  that  the  emerging  radiation  can  be  re- 
solved into  an  upwelling  branch  and  a downwelling  branch.  In  satellite  work, 
interest  lies  in  the  former  for  this  is  what  the  infrared  radiometer  on  board 
a satellite  presumably  "sees."  Then,  Eq.  (19)  may  be  decomposed  into  two 
equations. 


Upward: 


0 < p < 1 


d,I(ld7+~  ^ = I(t;  +u)  ' (1  ‘ “o}  Bv[T(t)] 


p(+p,  ip1 ) I(t;  ±p‘ ) dy1 


(21a) 


Downward: 


0 < p < 1 


dl(jdT"~~  = 'I(T;  “u)  + (1  • Bv[T(l)] 

to  /*  + l 

+ J-  I p(-y,  ±y')  I(t;  ±p ' ) dp1 

-1 


(21b) 


The  scattering  integrals  in  Eqs.  (21a,  b)  are  represented,  respectively,  by 


p(+P,  ±u') 


I(x;  ±p‘ ) dp' 


o 


p(+p,  +p' ) I (i ; + p1 ) 


p(+p , -p' ) I (t ; -p' ) 


dp' 


dp' 


(22a) 
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p(-y,  ±y' ) I(t ; ±y' ) dy 


p(-y,  +y‘ ) I(x ; +y 1 ) dy1 


0 


(22b) 


With  the  above  expressions,  the  derivation  of  the  radiative  transfer  equation 
for  the  infrared  spectral  region  is  complete.  The  solutions  of  Eqs.  (21a,  b) 
for  the  purpose  of  examining  the  optical  properties  of  clouds  in  satellite 
applications  will  be  presented  in  later  sections. 


EXTINCTION  COEFFICIENTS  AND  PHASE  FUNCTION 

In  the  radiative  transfer  equation,  there  are  several  optical  parameters 
that  need  be  fixed,  namely,  the  extinction,  scattering,  and  absorption 
coefficients,  B , B$,  and  6a,  derivable  from  the  so-called  Mie  parameters; 

the  optical  depth  t,  obtainable  directly  from  these  coefficients,  as  demon- 
strated in  Eq.  (5);  and  the  phase  function  which  is  usually  made  up  of  two 
scattering  functions,  Rayleigh  and  Mie.  In  infrared  radiation  where  thermal 
radiation  dominates  almost  completely  (see,  for  example,  Yamamoto  et  al . [6]), 
Rayleigh  scattering  by  air  molecules  is  generally  neglected,  although  it  can 
be  easily  computed.  It  may  be  included  in  the  following  manner: 


p(x;  y,  y')  = T(x)  M(y,  y')  + [(1  - T(t ) ] R(y,  y')  (23) 

where  T(t)  is  the  turbidity  factor  defined  by 


» 

1 

I 

1 

3 

i 

i 

l 

I 

i 

1 

| 


T(t) 


T /(t  + T ) 

s ,m  v s ,m  s ,r' 


(24) 


where  t and  t represent,  respectively,  the  optical  depths  due  to 
s jin  s » f 

Mie  and  Rayleigh  scattering.  In  Eq.  (23),  M(y,  y')  and  R(y,  y')  are 
the  Mie  and  Rayleigh  phase  functions,  respectively. 

However,  molecular  absorption  in  the  9-  to  12-micrometer  window  is  not  negli- 
gible. The  most  important  absorber  in  the  window  is  water  vapor.  As  a result, 
the  optical  depth  parameter  has  one  component  due  to  water  vapor  absorption 
in  addition  to  the  absorption  by  aerosol  and/or  cloud  particles. 
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Since  the  underlying  theory  and  the  computational  techniques  for  these 
parameters  including  atmospheric  transmission  have  been  discussed  in  some 
detail  by  Gomez  [7],  Gillespie  and  Petracca  [8],  and  Gomez  et  al . [9] 
among  others  in  published  ECOM  reports,  this  short  section  can  serve 
as  a footnote  to  them.  As  a matter  of  fact,  the  present  work  may  also 
be  looked  upon  as  a continuation  of  their  effort  in  target  signature 
studies . 

It  should  be  pointed  out  that  this  laboratory  has  several  computer  programs 
for  computing  the  Mie  coefficients  and  the  phase  function  for  any  given 
particle  size  distributions  and  complex  refractive  indices  and  for  calculat- 
ing high-resolution  atmospheric  transmission  due  to  molecular  absorption. 


CLOUD  MICROSTRUCTURE 

No  matter  how  elegant  mathematical  techniques  may  be  in  computing  the 
Mie  coefficients  and  in  constructing  the  Mie  phase  function  or  scatter- 
ing diagram,  results  from  the  solution  of  the  radiative  transfer  equation 
will  not  be  very  meaningful  if  the  particle  data  and  the  complex  refractive 
indices  used  in  computing  them  are  questionable.  Therefore,  it  is 
important  to  gain  a little  understanding  of  cloud  physics. 


Our  planetary  atmosphere  is  never  free  of  aerosol  particles,  but  it  is 
free  of  cloud  particles  for  only  a small  part  of  the  globe,  as  a glance 
over  the  GOES  visible  or  infrared  imagery  will  readily  show.  Aerosol  is 
a colloidal  system  in  which  the  dispersed  phase  is  composed  of  either  fine 
solid  or  liquid  particles  and  in  which  the  dispersion  medium  is  usually 
the  air.  There  is  no  clear-cut  upper  limit  to  the  size  of  particles 
comprising  the  dispersed  phase  in  an  aerosol,  but  as  in  all  other 
colloidal  systems,  it  is  rather  commonly  set  at  1 micrometer.  Haze, 
most  smokes,  and  some  fogs  and  clouds  may  be  regarded  as  aerosols. 

However,  it  is  not  good  usage  to  apply  the  term  to  ordinary  clouds  whose 
droplets  are  so  large  as  to  rule  out  the  usual  concept  of  colloidal 
stability.  It  is  poor  usage  to  apply  the  term  to  the  dispersed  particles 
alone,  and  aerosol  is  a system  of  dispersed  phase  and  dispersing  medium 
taken  together. 


With  the  aerosol  thus  clarified,  a cloud  can  now  be  defined.  A cloud 
is  a visible  aggregate  of  minute  water  and/or  ice  particles  in  the 
atmosphere  above  the  earth's  surface.  Cloud  differs  from  fog  only  in 
that  the  latter  is,  by  definition,  in  contact  with  the  earth's  surface. 
Clouds  form  in  the  free  atmosphere  as  a result  of  condensation  of 
water  vapor  in  rising  currents  of  air  and  hence  by  cooling  or  by  the 
evaporation  of  the  lowest  stratum  of  fog.  For  condensation  to  take 
place  at  the  point  of  saturation  or  a low  degree  of  supersaturation  in 
the  atmosphere,  condensation  nuclei  are  required  for  the  formation  of 
water  clouds,  while  ice  nuclei  are  required  for  ice-crystal  clouds.  The 
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size  of  cloud  droplets  varies  from  one  cloud  type  to  another.  Within 
any  given  cloud,  there  always  exists  a finite  range  of  sizes.  Generally, 
cloud  droplets  range  between  1 and  100  micrometers  in  diameter,  and  hence 
are  very  much  smaller  than  raindrops,  usually  of  the  order  of  1 millimeter 
or  greater. 

There  is  a connection  between  aerosol  and  cloud  particles,  but  the  connec- 
tion is  limited  only  to  those  aerosol  particles  which  are  in  the  neighbor- 
hood of  0.1  micrometer  or  larger  in  diameter  and  which  are  hygroscopic  in 
the  case  of  water  clouds  or  have  the  capacity  for  ice  nucleation  in  the 
case  of  ice  clouds.  These  aerosol  particles  which  have  the  ability  to 
serve  as  such  nuclei  are  generally  referred  to  as  cloud  condensation  nuclei 
or  ice  nuclei.  The  rest  of  the  aerosol  particles  that  do  not  partake  in 
the  cloud-  or  fog-forming  process  can  be  classified  as  haze  particles. 
Another  phenomenon  is  what  is  often  called  the  dry  haze,  the  visible  evi- 
dence of  the  existence  of  aerosol  in  the  air,  when  atmospheric  humidity  is 
low,  or  damp  haze  when  humidity  is  relatively  high.  In  the  latter  case, 
these  haze  particles  must  be  hygroscopic.  Near  or  at  saturation,  they  may 
then  form  mist. 

Given  a slight  supersaturation  of  the  order  of  perhaps  0.01%,  some  of 
these  haze  particles  or  cloud  embryos  may  overcome  the  critical  free 
energy  point  to  grow  to  larger  sizes.  In  a damp  haze,  there  could  be  a 
few  scattered  droplets,  but  they  are  of  insufficient  number  to  appear  as 
a cloud  or  fog.  Figure  3 is  known  as  the  Kohler  curve  of  droplet  growth. 

The  highest  point  in  the  curve  is  the  critical  point  which  designates 
the  so-called  critical  supersaturation  and  the  critical  radius  of  a 
condensation  nucleus.  Thus,  the  critical  point  is  the  dividing  line 
between  what  the  cloud  physicist  refers  to  as  the  haze  particles  on  the 
stable  side  and  the  cloud  droplets  on  the  unstable  side.  When  a condensa- 
tion nucleus,  necessarily  of  a sufficient  size  and  hygroscopici ty , grows 
into  a droplet,  it  may  be  said  that  the  droplet  has  reached  infinite 
dilution.  For  all  practical  purposes,  a cloud  droplet  can  be  considered 
a pure  water  droplet. 

The  somewhat  lengthy  description  of  aerosol,  haze,  and  cloud  particles 
given  above  was  taken  from  well-known  cloud  physics  books  [10,11,12]  and 
the  Glossary  of  Meteorology  [13]  in  order  to  clarify  that  aerosol,  haze, 
and  cloud  particles  do  not  mean  the  same  thing.  In  dealing  with  the 
aerosol  particles,  not  all  of  which  can  properly  be  regarded  as  spherical 
particles,  chemical  composition  should  be  of  concern.  In  contrast,  the 
haze  particles  and  the  cloud  droplets  can  be  treated  as  spherical  particles. 
In  the  former,  chemical  composition  is  still  important,  while  in  the 
latter  it  can  be  disregarded.  In  optical  property,  calculations  show 
that  the  extinction  coefficient  of  a spherical  particle  10  micrometers 
in  diameter  is  about  30  times  greater  and  that  of  a 5-micrometer  particle 
about  10  times  greater  than  that  of  a 1-micrometer  particle  in  the  10- 
micrometer  window  without  regard  to  their  composition.  When  particles 
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Figure  3.  The  equilibrium  supersaturation  as  a function  of  droplet 
radius  for  solution  droplets  containing  the  indicated 
masses  of  sodium  chloride.  (From  Mason,  1971) 
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range  from  0.0001  to  0.1  micrometer  in  diameter,  as  are  usually  found  in 
aerosol,  their  effect  on  transmission  in  the  infrared  window  is  quite  negli- 
gible in  comparison  with  the  effect  of  cloud,  haze,  or  even  dust  particles 
(in  a sandstorm),  as  can  be  readily  inferred  from  the  figure's  depicting 
extinction  and  hence  transmission  in  the  EC0M  report  by  Gomez  et  al . [9]. 

The  study  of  cloud  microstructure  is  only  one  facet  of  our  interest  in 
cloud  physics  which  embraces  not  only  the  investigation  of  the  condensation 
process  and  precipitation  physics  in  clouds,  but  also  radiative  transfer, 
optical  phenomena,  electrical  phenomena,  and  a wide  variety  of  hydro- 
dynamic  and  thermodynamic  processes  in  clouds.  In  this  section,  we  shall 
investigate  the  statistical  features  of  cloud  microphysics,  ignoring  its 
dynamic  processes  leading  to  the  formation  of  the  various  cloud  types  and 
the  techniques  for  sampling  cloud  droplets.  Cloud  and  fog  sampling  has 
been  in  effect  over  40  years.  Innumerable  droplet  data  have  been  collected. 
Yet,  knowledge  of  the  different  types  of  cloud  microstructure  is  still 
incomplete.  The  complexity  arises  from  the  fact  that  the  microstructure 
is  not  only  a function  of  cloud  type  but  also  of  geographical  location, 
position  within  the  cloud,  and  time  of  cloud  development.  Great 
effort  has  been  made  in  categorizing  these  clouds  by  means  of  different 
statistical  parameters;  however,  success  has  generally  been  limited.  Never- 
theless, it  may  be  of  some  interest  to  show  how  cloud  microstructure  varies 
with  these  factors.  Figure  4,  reproduced  from  Diem  [14],  shows  the  change 
as  a function  of  cloud  type.  Figure  5 from  Squires  [15]  shows  the  change 
with  geographical  location,  Fig.  6 from  Zaitsev  [16]  the  change  with 
height  level,  and  Fig.  7 from  Low  [17]  the  change  with  time. 

Of  interest  to  the  cloud  physicist  are  such  parameters  as  the  size 
range,  mean  radius,  mode  radius,  median  radius,  mean- volume  radius, 
and  liquid  water  content  (LWC)  which  is  expressed  in  unit  of  grams  per 
cubic  meter  (g  nr3).  Thus,  particle  size,  distribution,  concentration, 
thermodynamic  phase,  shape,  and  orientation  constitute  the  so-called 
microstructure  of  clouds.  Of  no  less  importance,  but  so  far  not 
easily  measurable,  is  the  percentage  of  supersaturation  attained  in  a 
cloud.  In  fact,  it  may  be  judged  to  be  the  most  important.  Given  the 
supersaturation  and  the  size  distribution  of  condensation  nuclei  together 
with  their  hygroscopic  property,  these  statistical  parameters  can  be 
readily  estimated  (e.g.,  [18]).  From  condensation  nucleus  measurement, 
Twomey  [19]  estimated  that  supersaturation  in  natural  clouds  may  be  as  low 
as  0.14%  and  as  high  as  1.1%,  depending  upon  the  ascent  rate  of  the  parcel 
of  saturated  air  and  the  number  of  nuclei  therein.  Following  the  same 
I approach,  but  employing  a measuring  technique,  Low  [20]  showed  that  the 

t supersaturation  attained  in  the  ordinary  radiation-advection  fogs  varied 

from  about  0.03%  to  about  0.10%,  mainly  due  to  radiation  cooling. 

Next,  the  liquid  water  content  shall  be  mentioned  briefly.  The  liquid 
water  content  can  be  calculated  from  the  number  concentration  of  cloud 
[ droplets  (if  sampled  correctly)  and  measured  with  an  instrument  such  as 
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Figure  4.  The  mean  droplet-size  distribution  of  various 
cloud  types.  (From  Diem,  1948) 
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5.  Histograms  of  the  percentages  of  droplet  samples  taken, 
not  all  at  the  same  geographical  location,  in  each  of 
five  cloud  types.  (From  Squires,  1958) 


19 


- 1 • 


. y».  -vk  - .1%-^  » ti 


3 


£ 0200h 


oioof- 


0000  L 


I — I — I — I — I I — . — > ^ • I 

0.10  0.15  0 20 

(1)  UQUEO  WATER  CONTENT  (g  cm'3) 

I . ■ . i . . - . I 

90  110  13  0 

(2)  MEAN  VOLUME  RADIUS  (>im) 


I 1 I I l i 

0 0 .01  02  03  04  05 

(3)  ESTIMATED  SUPERSATURATION  (%) 


Figure  7.  The  droplet-size  distribution,  liquid-water  content,  mean 
volume  radius,  and  supersaturation  as  functions  of  time. 
(From  Low,  1975) 
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the  hot-wire  LWC  meter.  In  clouds  it  can  vary  from  about  0.15  g m-3  in  a 
dry  continental  cumulus  cloud  [21]  to  as  high  as  3.9  g m-3  in  a highly 
convective  cumulus  congestus  [22].  The  other  cloud  types  apparently  fall 
somewhere  in  between  these  extremes.  As  to  the  fogs,  the  microphysical 
data  are  not  so  extensive.  Again,  the  fog  microstructure  depends  upon 
the  type  and  other  factors;  it  may  be  wet  or  dry.  Inland  fogs  may  have 
an  average  value  of  about  0.17  g m~3  (e.g.,  [23]),  and  coastal  fogs  an 
average  value  of  about  0.41  [20]. 

With  the  droplet  data  available,  an  analytical  distribution  or  density 
function  can  be  chosen  to  fit  these  data.  For  raindrop  size  distribution, 
the  Marshal  1 -Palmer  [24]  appears  to  be  generally  accepted.  For  cloud 
droplet  size  distribution,  a number  of  distribution  functions  have  been 
suggested,  notably.  Best's  exponential  distribution  [25],  the  lognormal 
distribution  [26],  the  gamma  distribution  [27],  and  the  Khrgi an-Mazi n 
distribution  [28]. 

Because  of  the  complex  nature  of  cloud  droplet  data,  these  data  can 

fit  a variety  of  standard  or  nonstandard  density  functions,  depending  upon 

where,  in  what  type  of  cloud,  at  what  level,  and  when  they  are  collected. 

It  would  perhaps  be  somewhat  pedantic  to  deliberately  choose  a more  compli- 
cated distribution  function  than  a simple  one.  In  an  exhaustive  comparison 
among  Best's,  the  lognormal,  and  their  own  distribution  functions  on  the 
basis  of  approximately  half  a million  stratocumul us  cloud  droplet  samples, 
Khrgian  and  Mazin  [29]  noted  that  the  lognormal  distribution  showed  a 
slightly  greater  accuracy  than  either  of  the  others.  Borovikov  [11]  felt 
that  the  complexity  of  working  with  the  lognormal  distribution  rendered 
it  difficult  to  handle  in  certain  cases  and  would  perhaps  more  than  offset 
the  advantage  gained  in  greater  accuracy.  Until  more  convincing  evidence 
dictates  otherwise,  the  author  of  the  present  paper  shall  choose  the 
Khrgian-Mazin  distribution  function  because  of  its  greater  simplicity.  As 
such,  the  function  is  given  by 

n(r)  = Ar2  exp(-br)  , (25) 

where  n(r)  is  the  number  of  droplets  per  cm3  in  size  r per  unit  radius 
interval,  and  A and  b are  constants.  This  distribution  may  be  seen  as  a 
modified  gamma  distribution  of  two  parameters,  a and  s,  with  a set  to  2 
and  6 = 1/b.  In  fact,  when  a is  set  to  0,  the  Marshall-Palmer  distribution 
for  raindrops  obtains. 

As  a probability  density  function,  the  constant  A in  Eq.  (25)  becomes 

A = b3/2  (26) 

and  the  other  statistical  parameters  calculated  from  the  moment  generating 
function  of  the  gamma  distribution  are  as  follows: 
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(27) 


where  N is  the  total  number  of  droplets  per  cm3,  and  rm  is  expressed  in 

micrometers.  Here,  the  density  of  a cloud  droplet  is  assumed  to  be  unity. 
However,  this  assumption  cannot  be  conveniently  made  in  the  case  of  haze 
particles,  as  already  noted. 


With  the  distribution  function  thus  fixed,  attention  is  directed  to  the 
types  of  clouds  to  be  dealt  with.  The  principal  types,  with  the  exception 
of  the  cirrus,  but  not  their  variations,  will  be  considered.  There  is  an 
extreme  paucity  of  cirrus  data,  and  no  less  serious  is  the  uncertainty 
about  the  scattering  properties  for  the  irregularly  shaped  ice  crystals 
of  the  cirrus.  In  fact,  any  clouds  whose  tops  reach  well  above  the  freezing 
level  may  contain  ice  crystals.  Single  scattering  properties  of  these 
nonspherical  ice  crystals  would  be  difficult  to  generate.  The  following 
table  summarizes  the  statistics  of  the  warm  clouds  which  will  be  used  in 
the  present  report: 


CLOUD  TYPES,  TEMPERATURES,  AND  STATISTICAL  PARAMETERS 


Type 

T 

(°C) 

r . 
min 

(ym) 

r 

m 

r 

max 

LWC  (g  m-3) 

N 

b 

NA 

St,  Sc,  and  Ac 

5 

1 

4-7 

20 

0.12-0.40 

125-200 

0.750 

42.18750 

Cumul us 
(continental ) 

1 

2 

5-9 

25 

0.20-0.50 

75-170 

0.600 

18.36000 

Cumul us 
(mari time) 

3 

2.5 

8-12 

30 

0.40-0.80 

50-85 

0.375 

2.24121 

Cumul us 

(fairly  strong) 

0 

3 

]0_-15 

40 

0.60-1.50 

25-65 

0.300 

0.87750 

Note:  Underlined  numbers  are  being  used  as  computer  inputs. 


Surface  Temperature  = 32°C 
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By  scanning  the  microphysical  data  documented  in  the  books,  particularly 
Borovikov's,  the  author  mentally  averaged  the  numbers  and  then  juggled  them 
so  as  to  come  up  with  what  he  believed  to  be  reasonable  ranges  for  the  middle 
to  high  clouds  in  midlatitude  summer  months.  The  cloud  temperatures  were 
obtained  in  the  same  manner  by  consulting  the  climatological  tables  (e.g., 
Upper-Air  Climatology  of  the  United  States  [30]).  With  these  cloud  parameters 
the  radiative  transfer  equation  can  be  solved  and  the  optical  properties 
and  emergent  radiances  in  the  infrared  window  region  can  be  determined. 


SOLUTION  OF  THE  RADIATIVE  TRANSFER  EQUATION 

There  are  a number  of  techniques  for  solving  the  radiative  transfer  equation, 
completely  or  approximately.  Among  them  are  the  so-called  exact  solution  for 
pure  Rayleigh  scattering  [31],  the  doubling  method  for  a homogeneous  medium 
[32,33,34],  the  iteration  method  [35,36,37],  the  discrete  ordinate  method  [38], 
the  spherical  harmonics  technique  (e.g.,  [39]),  the  invariant  imbedding 
method  [40],  and  the  Monte  Carlo  method  [41].  These  different  methods  have 
been  reviewed  by  Yamamoto  and  Tanaka  [42]  and  by  Hansen  and  Travis  [43]. 

A majority  of  the  authors  dealing  with  transfer  problems  are  either  astro- 
physicists or  physicists.  Generally,  they  are  concerned  with  the  stellar  or 
solar  radiation.  Apparently,  they  were  not  interested  in  radiative  transfer 
in  the  infrared  spectral  region  until  the  successful  orbiting  of  the  TIROS 
series  in  early  1960  made  infrared  cloud  imagery  available  to  the  meteoro- 
logists day  and  night.  Even  so,  comparatively  fewer  papers  have  been 
published  on  this  subject.  Yamamoto  et  al.  [6]  seem  to  be  among  the  earliest 
atmospheric  physicists  to  deal  seriously  with  the  scattering,  absorptive,  and 
emissive  properties  of  water  clouds  in  the  10-micrometer  infrared  window, 
first  by  expanding  the  transmission,  scattering,  and  emission  functions  into 
power  series  and  then  by  iteration  to  obtain  the  final  solution. 

In  a later  paper,  Yamamoto  et  al . [44]  showed  how  to  obtain  integrated 
upwelling  and  downwelling  radiation  over  a broad  infrared  band  on  the  basis 
of  the  same  scheme.  Zdunkowski  and  Choronenko  [45]  employed  a straightforward 
iterative  technique  to  investigate  the  blackness  of  isolated  clouds.  Hunt 
[46]  used  the  matrix  operation  method  to  solve  the  transfer  equation  in  both 
visible  and  infrared  spectral  regions  and  presented  numerous  tables  to  depict 
the  optical  properties  of  both  water  and  ice  clouds;  in  doing  so,  he  assumed 
the  ice  crystals  to  be  spherical.  In  all  these  studies,  clouds  are  assumed 
to  be  homogeneous  in  both  temperature  distribution  and  droplet  number  concentra- 
ti  on. 

Below,  the  equation  will  be  solved  by  the  Gauss-Seidel  iteration  method  [47], 
following  Herman  and  Browning  [35].  Equations  (21a)  and  (21b)  are  reformulated 
as  follows: 


where 


I(t;  +p)  - 0 ( t ; +p) 


y Ui^dx  ~ ^ “ ’ I(t;  -y)  + J(t;  -y) 


J(r;  +p)  = (1  - W ) B[T(t ) ] +2 - I p(+p;  ±p')  I ( T ; ±p')dp‘ 

-1 


J(t;  -p)  - (1  - wQ)  B[T(  t ) ] + y~  I p(-p;  ±p')  I(t;  ±p')du' 


Note  that  the  scattering  integrals  are  given  by  Eqs.  (22a)  and  (22b). 

Figure  1 shows  that  the  cloud  may  be  sliced  into  N thin  layers  of  equal 
geometric  thickness,  as  in  our  case,  or  of  equal  optical  thickness,  from  the 
top  to  the  bottom  in  opposite  direction  to  geometric  height;  thus,  the  upper 
boundary  of  the  cloud  is  labeled  n = 1,  and  the  first  layer  n = 1 for  n = 1, 
2,  3,  . . . N + 1 levels  and  n = 1,  2,  3,  . . . N layers.  For  the  present, 
the  optical  mass  from  space  down  to  the  upper  boundary  (i.e.,  n = 1)  is  taken 
to  be  nil;  hence,  the  optical  thickness  at  n = 1,  =0.  Then,  T2  = Ti  + 

t2,/t3  = ti  + t2  + x3’  and  so  on  ^or  ^eve1s  2,  ...  N + 1.  With  this  in 

mind,  multiply  Eq.  (21a)  by  the  integrating  factor,  e+t^'  and  Eq.  (21b)  by 

the  factor,  e-t^u,  and  then  integrate  the  resulting  expressions  successively 
from  x -|  to  t3,  t 2 to  and  so  on.  The  following  iteration  formulas  are 

obtai ned: 


I(xn;  +p)  = I(Tn+1;  +p)  e-(jn+rTn)/p  + f " ^(t;  +u)  e'(t_Tn)/v  y 


(28a) 


I ( t n ; -y)  = I(xn_.| ; -p)  e"^Tn"Tn-l  + J " J(t;  -p)  e"^Tn_t^v 


(28b) 


For  relatively  thin  layers  of  the  order  of  perhaps  20  to  50  meters  thick, 
depending  on  the  type  of  cloud  studied,  an  arithmetic  mean  of  the  source 
function  J may  be  taken.  Now  let  J(xn;  ±y)  be  the  average  source  function 

at  the  top  and  bottom  of  the  nth  layer  and  At.  = t - t , where  k varies 

k n n-1 

with  the  specification  of  n.  Equations  (28a)  and  (28b)  may  then  be  put 
into  a rather  compact  form  as  follows,  upon  integration: 

I (Tn » tv)  = 1 (Tn±i  ’ tv)  e AT|</^  + T(xk;  ±u)  (1  - e'ATk/u)  , (29) 

where  positive  values  of  y (i.e.,  +y,  upwardgoing) 

k = n , and  1 _<  n < N + 1 , 

and  for  negative  values  of  y (i.e.,  -y,  downwardgoing) 

k = n - 1 , and  1 < n <_  N + 1 

In  using  these  compact  notations,  however,  it  should  be  noted  that  J ( xn ; +y) 
and  J(xn;  -y)  are  taken  to  mean  the  average  radiation  emerging  at  the  top 
of  the  nth  layer  and  at  the  bottom  of  the  (n  - l)th  layer,  respectively. 
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At  this  point,  it  may  be  appropriate  to  introduce  the  boundary  conditions, 
noting  that  radiative  transfer  is  a two-point  boundary  value  problem.  It 
is  the  usual  practice  in  the  literature  to  assume,  though  not  quite  correctly, 
that  in  the  infrared  region  there  is  no  radiation  incident  on  the  top  of 
the  cloud.  This  assumption  shall  be  adopted  for  the  present;  then 


I ( t ^ ; -\i)  = 0 , 


(31a) 


I 
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where  refers  to  the  optical  thickness  from  space  to  the  top  of  the 
cloud,  which  has  already  been  assumed  to  be  zero. 

At  the  bottom  of  the  cloud,  if  an  isolated  cloud  is  investigated  as  in 
the  paper  by  Zdunkowski  and  Choronenko  [45],  the  following  would  apply 


I(tb;  +u)  = 0 (31b) 

On  the  other  hand,  the  ground  surface  may  be  the  boundary  in  the  case 
of  a dry  column  below  the  cloud  base;  or  a humid  column  together  with 
ground  radiation  below  the  cloud  base  may  be  the  boundary.  In  this 
present  study,  the  former  is  adopted;  thus. 


I (Tb ; +p)  = B(Tg)  , (31c) 

where  T is  the  temperature  of  the  ground  surface  whose  emissivity  is 
assumed9to  be  1.  Since  this  discussion  is  more  centered  in  the  optical 
properties  of  various  cloud  types  at  present,  such  as  cloud  emissivity 
and  transmissivity,  Eqs.  (31a),  (31b),  and  (31c)  will  be  used. 

Either  the  top  level  (i.e.,  n = 1)  or  the  bottom  level  (i.e.,  n = N + 1) 
can  be  used  as  a start.  For  this  discussion,  begin  with  the  top  level  and 
proceed  to  the  base.  During  the  downward  excursion  through  each  successive 
layer  to  the  base,  assume  that  there  shall  be  no  upward  radiation.  However, 
in  going  from  the  first  to  the  second  layer,  the  scattering  due  to  downward 
emission,  given  by  the  second  term  on  the  right  side  of  Eq.  (22b)  can  no 
longer  be  neglected;  it  must  be  included  in  the  computation.  As  the  steps 
are  traced  upward  through  each  succeeding  layer,  Eqs.  (22a)  and  (22b)  come 
into  the  picture.  When  the  top  level  is  reached,  the  first  iteration  has 
been  completed.  Equivalently,  it  may  be  said  that  the  upward  and  downward 
radiation  from  first  order  scattering  has  been  obtained.  As  the  computation 
undergoes  each  successive  iteration,  each  successively  higher  order  scat- 
tering is  taken  into  account.  The  iteration  procedure  terminates  when  for 
m > l 
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(32) 


!m+l  ^tN+1  ; ; ’ e ’ 

where  e has  been  set  to  1.002  in  this  study. 

So  far  the  computation  of  the  indispensable  phase  function  has  not  been 
mentioned.  Two  separate  computer  programs  have  been  obtained  through  the 
courtesy  of  NASA  for  the  computation  of  the  function.  One  is  slightly 
more  efficient  than  the  other;  however,  the  radiative  transfer  computer 
program  used  will  accept  the  coefficients  of  the  Legendre  polynomials 
given  by  any  phase  function  computations.  Dei rmendj i an 's  book  [48]  is 
useful  for  a description  of  the  computational  technique  for  constructing 
the  Mie  scattering  functions. 

In  concluding  this  section,  it  may  be  pointed  out  that  the  iteration 
scheme  as  outlined  above  will  render  it  possible  to  calculate  the 
transfer  of  infrared  radiation  in  nonhomogeneous  atmospheres  or  clouds. 

It  is  also  possible  to  consider  the  variation  of  the  distribution 
functions  of  molecules  and/or  particulates,  although  a large  amount  of 
computer  time  may  be  required. 


NUMERICAL  RESULTS 

Numerical  solutions  of  complex  integro-differential  equations  invariably 
present  the  question  of  how  reliable  and  accurate  the  numerical  results  are 
despite  careful  attention  to  the  formulation  of  the  algorithm  and  the  con- 
vergence criterion.  The  best  way  to  resolve  the  question  is  to  compare  the 
results  with  controlled  or  selected  experiments.  The  authors  of  the  present 
project  hope  to  conduct  such  an  experiment  in  the  near  future.  A number  of 
radiation  results  generated  by  various  numerical  techniques  have  already 
been  published.  Thus,  from  the  comparisons  of  these  methods,  definitive  in- 
formation may  be  derived  concerning  the  reliability  of  each  method  employed 
in  the  calculations. 

A comparison  of  the  intensity  values  given  in  graphical  form  by  Yamamoto 
et  al . [6]  was  more  tractable  because  of  their  utilization  of  Dei rmendj ian ' s 
cumulus  cloud  model  [49]  whose  phase  functions  were  listed  for  several  wave- 
lengths in  the  infrared  band.  Therefore,  it  was  possible  first  to  check  the 
phase  function  of  this  project  for  the  same  cloud  model  against  that  of 
Dei rmendj i an' s (as  later  tabulated  in  his  book  [48])  and  then  to  compare 
intensity  values  with  Yamamoto's  at  the  cloud  top  and  other  levels,  utiliz- 
ing their  model  cases.  The  intensity  values  carefully  extracted  from  their 
graphs  were  found  to  lie  well  within  5%  of  ours  in  zenith  angles  from  0 up 
to  about  60  degrees.  Comparison  with  Hunt's  cloud  emissivity  and  trans- 
missivity values  [46],  however,  shows  some  differences  because  of  the  uncer- 
tainty involved  in  the  values  of  phase  functions.  Nevertheless,  for  the 
two  cases  of  water  cloud  at  the  11-micrometer  wavelength,  emissivity  and 
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transmissivity  values  of  this  project  lie  within  10%  of  his  tabulated  values 
using  Hunt's  cloud  models.  Therefore,  it  may  be  stated  with  confidence  that 
the  numerical  outputs  from  this  computer  program  employing  the  iteration 
method  are  at  least  comparable  with  those  given  by  the  adding  method  [46] 
and  the  technique  of  scattering  and  transmission  functions  [6]. 

As  input  data  to  the  computer  program,  the  values  of  the  various  cloud  types 
presented  in  the  table  were  used  and  the  complex  refractive  indices  for  water 
sphere  compiled  by  Irvine  and  Pollack  [50]  in  the  9-  to  12.5-micrometer  band, 
which  approximately  corresponds  to  the  bandwidth  of  infrared  imagery  in  cur- 
rent use.  Figure  8 shows  the  size  distributions  of  three  cloud  models.  The 
droplet  size  and  the  concentration  scales  for  the  stratus  are  at  the  top  and 
on  the  right  side  of  the  figure,  respectively.  The  stratus  model  is  again 
used  to  demonstrate  in  Figure  9 how  its  volume  extinction  coefficient  and  the 
single  scattering  albedo  vary  with  wavelength.  This  illustration  implies 
that  optically  a cloud  behaves  differently  at  different  frequencies.  The 
same  figure  also  shows  the  change  of  blackbody  radiation  with  frequency.  The 
optical  properties  of  other  cloud  types  will  be  different,  but  their  varia- 
tions with  wavelength  follow  closely  the  curves  as  depicted  in  Figure  9. 

They  will  not  be  repeated  for  each  case.  Figure  10  is  presented  as  an  example 
of  the  change  of  cloud  emissivity  (or  transmissivity)  with  wavelength  and 
cloud  depth.  The  cloud  flux  emissivity  is  usually  defined  by 


eA(t) 


I x ( r ; p ) pdu/Bx(Tc) 


(33) 


whereas  from  a satellite  point  of  view  the  upwelling  cloud  emissivity 
may  be  expressed  by 


- yo'-'i'MV 


(34) 


where  T is  the  cloud  temperature  B (T  ) its  blackbody  intensity,  and 

I (0,1)  the  upwelling  radiance  at  the  cloud  top.  A similar  expression 

holds  for  cloud  transmissivity,  except  that  the  cloud  temperature  is 
now  replaced  by  the  Planck  function  of  the  surface  temperature  and  the 
assumption  is  made  that  there  is  no  cloud  emission. 


Since  the  present  interest  is  in  infrared  imagery,  only  the  11.5-micrometer 
window  will  be  investigated,  to  conserve  computer  time.  The  11.5-micrometer 
window  lies  approximately  in  the  middle  of  the  frequency  response  curve 
of  infrared  imagery.  Figure  11  shows  how  the  emissivity  of  each  cloud 
model  varies  with  cloud  depth.  Ordinarily,  infrared  imagery  at  or  near 
the  subsatellite  point  is  of  more  importance  in  the  determination  of 
cloud-top  height  since  such  imagery  enables  the  determination  to  be  made 
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with  greater  confidence.  For  this  reason.  Figure  12  depicts  the  normal 
upwelling  radiation  at  the  zenith,  as  presumably  observed  by  a satellite 
radiometer.  The  insert  shows  an  enlargement  of  the  segment  where  the 
stratus  emission  approaches  blackbody  emission  at  about  650  meters  above 
the  cloud  base  when  no  emission  from  the  surface  is  considered.  When  the 
surface  emission  is  included,  the  cloud  appears  to  radiate  more  energy 
at  650  meters,  thereby  causing  it  to  become  "blacker"  at  that  level  than 
the  blackbody  at  the  cloud  temperature.  The  normal  upwardgoing  intensities 
at  the  zenith  of  other  cloud  types  are  presented  in  Figure  13. 


DISCUSSION 

In  the  present  investigation,  four  cloud  models  have  been  chosen  which 
are  believed  to  be  representative  of  the  four  cloud  types.  After  experi- 
mentation with  the  stratus  cloud  of  various  thicknesses,  it  was  discovered 
that,  given  the  microphysical  parameters  of  the  cloud  as  tabulated  in  the 
table,  no  thickness  greater  than  1 kilometer  needed  to  be  considered. 
Moreover,  the  level  at  which  the  same  cloud  approaches  the  property  of  a 
blackbody  varies  with  wavelengths,  as  can  be  readily  inferred  from  the 
variation  with  wavelength  of  the  volume  extinction  coefficient  in  Figure  9. 
The  same  cloud  will  appear  the  "blackest"  in  the  11-  to  13-micrometer  window 
region,  which  implies  that  the  thickness  required  for  it  to  radiate  like  a 
blackbody  will  be  at  the  minimum.  In  fact,  as  shown  by  Hunt  [46],  for  t = 
100,  the  emissivity  of  the  same  cloud  is  0.8980  in  the  3.8-micrometer  window, 
but  0.9984  in  the  1 1 -mi crometer  window.  To  examine  the  internal  radiation 
in  detail,  the  cloud  has  been  divided  into  50  sublayers.  However,  this  is 
not  necessary  in  practice  if  interest  is  only  in  the  upwelling  radiation 
emerging  from  the  cloud  top.  The  test  runs  reveal  that  upwelling  radiances 
vary  within  0.2%  by  employing  either  10  sublayers  of  100  m each  or  50  of  20 
m each  for  a cloud  layer. 

Figure  10  again  indicates  that  the  stratus  approaches  blackbody  radiation 
at  a lower  level  at  12.5  than  at  9 micrometers.  However,  at  the  same  11.5- 
micrometer  wavelength,  as  depicted  in  Figure  11,  the  cloud  with  greater 
liquid  water  content  attains  blackbody  radiation  at  a lower  level  than 
the  one  with  smaller  liquid  water  content.  Thus,  while  it  takes  the 
stratus  about  550  m in  thickness  to  achieve  blackbody  radiation,  a fairly 
strong  convective  cloud  requires  only  about  220  m,  about  half  as  much. 

, In  terms  of  liquid  water  content,  the  latter  cloud  contains  five  times 

« as  much  water  as  the  former.  Because  of  the  tremendous  amount  of  computer 

time  required,  it  has  not  been  possible  to  establish  a tractable  relation- 
ship between  the  liquid  water  content  of  a cloud  and  the  thickness  required 
i of  the  cloud  for  blackbody  radiation.  Returning  to  Figure  10,  it  may  be 

■ noted  that  at  either  9 or  12.5  micrometers  the  sums  of  flux  emissivity  and 

» transmissivity  at  all  depths  are  close  to  unity;  this  implies  that  in  the 

i 10-micrometer  infrared  band  there  is  little  or  no  flux  reflectivity  of 

i concern  since 
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NORMAL  UPWARD  INTENSITY  (mW  m sr  cm)  FROM  200  TO  500  METERS 


t + e + r = 1 


(35) 


where  r is  the  reflectivity  by  the  warm  clouds. 

When  a cloud  reaches  a thickness  such  that  it  radiates  almost  like  a 
blackbody,  the  upwelling  radiation  coming  from  the  cloud  top  becomes  very 
nearly  isotropic,  except  towards  the  limbs.  As  was  also  found  in  our 
calculations,  Zdunkowski  and  Choronenko  [45]  showed  isotropic  radiation 
at  the  cloud  top  for  the  case  of  an  isolated  cloud  (i.e.,  without  consider- 
ing surface  emission)  and  Yamamoto  et  al . [6]  for  the  cases  of  nonisolated 
clouds  (i.e.,  with  surface  radiation  considered).  Since  infrared  imagery 
is  gray-scale  contours  of  cloud-top  intensities  (i.e.,  radiances),  the 
emergent  radiation  at  or  near  the  zenith  is  of  particular  interest,  as 
already  noted.  In  Figures  12  and  13,  attention  is  directed  toward  the 
intersections  of  the  solid  and  dashed  lines.  Although  only  the  upwelling 
radiation  has  been  displayed  at  the  zenith  for  each  case,  it  may  be 
mentioned  that  up  to  60  to  70  degrees  the  intersections  are  located  at 
the  same  level  for  the  same  cloud  model  as  a result  of  being  isotropic. 

When  there  is  no  surface  radiation,  the  solid  curve  approaches  a straight 
line  (i.e.,  blackbody  radiation)  at  a level  somewhat  lower  than  that  of 
the  intersection.  This  is  more  prominent  in  the  stratus  model.  As  the 
liquid  water  content  or  the  number  of  larger  droplets  increases,  the  effect 
of  surface  emission  diminishes.  The  isolated  stratus  reaches  blackbody 
radiation  at  about  660  m from  its  base,  a continental  cumulus  at  about 
450  m,  a maritime  cumulus  at  about  370  m,  and  a fairly  strong  cumulus  at 
about  330  m.  In  contrast,  when  surface  radiation  is  present,  the  inter- 
sections, beyond  which  the  cloud  will  radiate  as  a blackbody  at  its  own 
temperature,  are  located,  respectively,  at  about  940  m,  500  m,  390  m,  and 
330  m.  Only  in  the  last  case,  the  solid  and  dashed  lines  meet  at  nearly 
the  same  level  or  depth,  indicating  that  surface  radiation  has  little 
effect  on  the  level  at  which  a heavily  water-laden  cloud  attains  blackbody 
radiation,  whether  in  isolation  or  not. 

What  has  just  been  discussed  appears  to  be  in  conflict  with  the  results 
described  in  Figure  11  with  respect  to  ^oud  emissivity.  It  should  be 
pointed  out  here  that  this  inconsistency  depends  upon  the  application  of 
the  concept  of  emissivity.  The  conventional  concept  of  emissivity  is 
defined  on  the  basis  of  radiation  flux,  expressed  by  Eq.  (33).  But  from 
the  point  of  view  of  satellite  radiometer  measures,  the  emissivity  value 
should  be  given  by  Eq.  (34).  Thus,  it  is  important  that  cloud  emissivity 
be  examined  carefully  in  conjunction  with  satellite  applications.  The  cloud 
cases  investigated  here  show  that  the  stratus  reaches  a depth  of  about  940  m 
(in  theory)  instead  of  540  m before  it  becomes  a blackbody  from  radiometric 
requirements.  The  continental  cumulus  must  attain  about  500  m instead  of 
280  m,  the  maritime  cumulus  about  390  m instead  of  240  m,  and  the  fairly 
strong  convective  cloud  about  330  m instead  of  220  m. 
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Of  course,  no  cloud  is  ever  isolated,  nor  is  any  cloud,  except  the  fogs, 
ever  in  immediate  contact  with  the  ground  surface  without  some  intervening 
water  vapor  and  temperature  distributions,  as  are  assumed  in  our  present 
investigation  as  well  as  in  the  studies  by  Yamamoto  et  al . [6],  Zdunkowski 
and  Choronenko  [45],  and  Hunt  [46].  Nevertheless,  the  calculations  presented 
here  concerning  the  gross  optical  properties  of  clouds  are  by  no  means 
invalid,  but  rather  represent  a first  step  toward  a better  understanding 
and  a more  detailed  examination  of  radiative  transfer  in  clouds  as  is 
applicable  to  satellite  problems. 


CONCLUSIONS  AND  RECOMMENDATIONS 

Thicknesses  or  depths  required  for  blackbody  radiation  have  been  specified 
for  the  cloud  types  studied.  Below  these  depths,  the  clouds  will  appear 
warmer  than  they  actually  are,  as  evident  in  Figure  13.  When  the  cloud 
temperature  is  known,  it  appears  possible  to  deduce  the  surface  temperature 
in  the  absence  of  water  vapor  below  the  cloud  base. 

For  satellite  applications  in  which  upwelling  radiance  (i.e.,  intensity) 
is  measured,  the  conventional  concept  of  emissivity  can  no  longer  be 
applicable.  Instead,  emissivity  calculated  on  the  basis  of  upwelling 
radiance  should  be  introduced,  thus  making  it  appropriate  for  a proper 
interpretation  of  cloud-top  temperatures.  Note,  however,  the  greater 
the  amount  of  liquid  water  is  in  a cloud,  the  less  is  the  difference 
between  the  two. 

Reflectivity  by  warm  clouds  in  the  10-micrometer  infrared  band  is  normally 
negl igible . 

The  optical  properties  of  a cloud  are  functions  of  wavelength,  as  also  noted 
by  other  authors.  The  cloud  appears  more  "black"  in  the  neighborhood  of 
the  11-  to  13-micrometer  window  than  in  any  other  infrared  windows. 

The  liquid  water  content  of  a cloud  is  an  important  microphysical  parameter 
in  relating  directly  to  the  degree  of  blackbody  radiation.  It  requires 
further  investigation  to  establish  a functional  relationship  between 
the  liquid  water  content  of  the  cloud  and  the  associated  thickness  for 
blackbody  radiation. 

During  the  next  project,  the  computer  program,  which  is  now  capable  of 
solving  the  radiative  transfer  equation  for  an  inhomogeneous  cloudy 
atmosphere  in  the  absence  of  water  vapor  absorption,  will  be  modified  to 
include  such  absorption  both  below  and  above  the  cloud  as  well  as  the 
frequency  response  function  of  the  satellite  radiometer  in  current  use. 

In  so  doing,  it  will  be  possible  to  compare  theoretical  radiance  values 
with  those  observed  values  from  satellite  data  streams  when  the  radio- 
meter's field  of  view  is  partly  or  completely  covered  by  clouds. 
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